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Abstract: Sequences from a region of the nuclear ribosomal 18S rDNA gene of approximately 1 912 base pairs (bp) 
were used to generate a molecular phylogeny for the Pentatomomorpha based on 53 species representing 21 putative fami- 
lies. Phylogenetic analyses using the most parsimony method (MP), maximum likelihood method (ML), and neighbor 
joining method (NJ) showed strong support that the Pentatomomorpha lineage is a monophyly and the superfamily 
Aradoidea is a sister group to the remainder of the Pentatomomorpha (Trichophora). The Trichophora could be divided into 
two clades: one clade consisted of the monophyletic superfamilies Pentatomoidea and Pyrrhocoroidea; the other was mainly 
the polyphyletic superfamilies Lygaeoidea, Coreoidea and Idiostoloidea. The superfamilies Lygaeoidea and Coreoidea were 
both polyphyletic. Within Lygaeoidea, Piesmatidae was sister to Berytidae. They formed a clade locating at the basal of 
the Trichophora and distantly related to the other two families Lygaeidae and Rhyparochromidae. This research suggested 
that 185 rDNA was a proper marker to reconstruct the phylogeny of Pentatomomorpha that was accordant to morphological 
studies and the research of Li et al (2005). The Pyrrhocoroidea was further divided from the Coreoidea (s./at). It was 
suggested that the Piesmatidae might be assigned as a superfamily of Pentatomomorpha rather than a family in Lygaeoidea. 
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The Pentatomomorpha is the second largest subor- der and one of the most important groups of the Het- 
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eroptera, numbering from 12 500 to 15 000 species 
worldwide. Of the 29 families of Heteroptera (Schuh & 
Slater, 1995) , 20 have been recorded in China. All of 
the taxa in this infraorder are terrestrial, and most are 
plant-feeding. They are important in agriculture, as 
some of them are pests. The Pentatomomorpha has an 
important phylogenetic position because it is one of the 
most phylogenetically distal groups in the Heteroptera. 
Since the twentieth century, many comparative mor- 
phology studies have provided strong evidence and im- 
portant information for the construction of a reasonable 
classification system for the Heteroptera. Tullgren 
(1918) established Trichophora based on the abdomi- 
nal trichobotheia, which was the origin of * Pentatomo- 
morha” in the present classification system. Leston et 
al (1954) established the infraorder Pentatomomorpha 
based on the pretarsal structure, male and female geni- 
talia; wing venation, egg morphology, and salivary 
gland structures. The superfamily Aradoidea (including 
Araridae and Termitaphididae ) and the group Tri- 
chophora make up the present Pentatomomorpha (Hen- 
ry, 1997). Schuh (1979) and Wheeler et al (1993) 
had proved that Pentatomomorpha was a monophyly us- 
ing cladistic analysis, and the sister group appeared to 
be Cimicomorpha. 

Up to 
(Aradoidea + (Pentatomoidea + the remainder of Tri- 


now, the phylogenetic hypothesis 
chophora) ) of this infraorder has been supported by Le- 
ston (1958) and Henry (1997), and accepted by most 
researchers. Within the Pentatomomorpha, the relation- 


Only the 


Aradoidea and Pentatomoidea are consistently recog- 


ships of superfamilies are uncertain. 
nized as monophyletic within these schemes (Henry, 
1997). Studies have recognized four ( Schaefer, 
1964), five (Sys, 1961, 1967; Schaefer, 1993; 
Schuh & Slater, 1995) , six (Carver et al, 1991; Hen- 
ry, 1997), or seven (Schuh, 1986; Henry & 
Froeschner, 1988; Li & Zheng, 1994) superfamilies. 
No consensus exists for the relative position of the Pies- 
matidae and Idiostolidae. The other superfamilies 
Coreoidea, Lygaeoidea, and Pyrrhocoroidea have been 
defined with some difference. More recently, Henry 
(1997) brought together all the variable morphological 
character information found in the literature to analyze 
it cladistically, and recognized six superfamilies 
( Aradoidea, 


coroidea, Idiostoloidea, Lygaeoidea ). The Piesma- 


Pentatomoidea, Coreoidea, Pyrrho- 
toidea was placed into the Lygaeoidea. Until now, it 
was showed that the relationships among the superfami- 


lies in Pentatomomorpha and the families within the su- 


perfamilies were still uncertain. 

With further research into insect molecular phylo- 
genetics, great progress has been made using molecular 
markers. The genes of rDNA and mtDNA are the two 
most widely used markers. Thus far, little attempt has 
been made to incorporate molecular data into phyloge- 
netic analyses of the Pentatomomorpha. Xie et al 
(2005) used 18S rDNA to recover phylogeny for the 
main lineages 


( Pentatomoidea, Coreoidea, Pyrrho- 


coroidea, and Lygaeoidea ) of Trichophora with 
Bayesian analysis, using Pentatomoidea as an out- 
group. The results showed that the relationship of 
Coreoidea, Pyrrhocoroidea, and Lygaeoidea was 
Pyrrhocoroidea + ( Coreoidea + Lygaeoidea), which 
differed to the hypothesis of Henry (1997) and Leston 
(1958) that was (Coreoidea + Pyrrhocoroidea) . Li et al 
(2005) were the first to analyze the whole Pentatomo- 
morpha phylogeny using 18S rDNA and the mitochon- 
drial COX1 gene based on 40 species representing 17 
putative families. The phylogenies were mostly congru- 
ent with morphological studies. Results strongly sup- 
ported the monophyly of Pentatomomorpha, and the 
placement of Aradoidea as sister to Trichophora. The 
monophyletic Trichophora was grouped into two major 
lineages: one being the monophyletic Pentatomoidea , 
and the other comprising Lygaeoidea, Coreoidea, and 
Pyrrhocoroidea. The phylogeny of the whole Pentatomo- 
morpha was: Aradoidea + (Pentatomoidea + (( Pyrrho- 
coroidea + Coreoidea + Lygaeoidea) ) ). Thus the results 
supported that the morphological characteristic abdomi- 
nal trichobothria was a synapomorphy of Trichophora. 
This is significant to determining the polarities of mor- 
phological characteristics in phylogenetic studies of 
Pentatomomorpha. Li et al (2005) suggested that 18S 
rDNA may be an optimal molecular marker for the phy- 
logeny of Pentatomomorpha; on the contrary, the 
COXI segment might not be. The position of Piesmati- 
dae, and the relationships among Pyrrhocoroidea, 
Coreoidea, and Lygaeoidea will require further sam- 
pling and more complete studies to resolve. 

The 18S rDNA gene is one of the most widely 
used markers in studies of insect phylogeny. Nowa- 
days, many sequences of Pentatomomorpha 185 rDNA 
have been accumulated in GenBank. According to the 
research of Li et al (2005), we chose 18S rDNA gene 
to generate the phylogeny of Pentatomomorpha. The 
study included additional superfamilies and families, 
which were Trisecus pictus (Idiostoloidea) , Canopidae, 
Cydnidae, and Hyocephalidae, based on 53 species 
representing 21 putative families and six superfamilies . 
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( MP), 
(ML), and neighbor joining (NJ) methods were used 


The most parsimony maximum likelihood 
to reconstruct the relationships among superfamilies in 
the Pentatomomorpha lineage. The position of Piesmati- 


dae is also discussed. 


1 Materials and Methods 
1.1 Taxa 


Fifty three taxa representing six superfamilies and 
2] putative families were analyzed to study the relation- 
ship at the high level of Pentatomomorpha. The speci- 
mens used and their current taxonomy are shown in 
Tab. 1. 

We collected 25 species in Guangdong Province 
during 2003 and sequenced their 18S rDNA. The nu- 
cleotide sequence data has been deposited in GenBank 
and their accession numbers are listed in Tab. 1. Other 
sequences were collected directly from GenBank. Three 
other species ( Emesaya brevipennis, Lygus hesperus 
and Campyloneura virgula) were used as outgroups be- 
longing to Cimicomorpha. Accession numbers of these 
sequences are also listed in Tab. 1. 

1.2 Total genomic DNA isolation 

Total genomic DNA was extracted from frozen, 
ethanol-preserved (95% ethanol) specimens and from 
pinned-dried specimens, usually obtained from thoracic 
muscles. Heads, abdomens and legs were stored in 
70% alcohol as voucher specimens. The method of 
DNA extraction is referred to Wen & He (2003), and 
was revised. The protocol was as follows: each speci- 
men was softly ground in 100 pL buffer (10 mmol/L 
Tris, 1 mmol/L EDTA, 0.1 mol/L NaCl, pH 8.0). 
After adding Proteinase K (200 pg /mL), the ho- 
mogenate was incubated at 56 °C for 1 — 2 hours, 95 °C 
for 45 seconds, and then centrifuged for 5 minutes. 
The supernatant was collected as the template for the 
PCR. DNA from some specimens, especially a few 
pinned-dried specimens, was extracted with EaZy Nu- 
cleic Acid Isolation Kit (Omega Bio-tek) . 

1.3 18S rDNA sequence cloning and sequencing 

The primers are listed in Tab. 2. Primers were 
derived from Loxdale & Lushai (1998). PCRs were 
conducted in 25 pL volume containing 2 - 6 pL of 
DNA, 0.5 U Taq polymerase (Takara, Dalian), 1.5 
uL 10 pmol/L primers, 2.0 uL 2. 5 mmol/L dNTPs 
(Takara), and 2.5 uL 10 x buffer. Amplification con- 
ditions were 95 *C for 5 min; 35 cycles of 95 *C for 30 
s, 50 C for 40 s, 72 C for 2 min; 72 °C for 8 min. 
Target products from the PCR were purified with A- 
garose Gel DNA Purification Kit (Takara) and se- 


quenced directly. DNA sequencing was performed in a 
PE/ABI 377 automated sequencer. 
1.4 Sequence alignment and phylogenetic analyses 

The sequences were aligned using default parame- 
ters of Clustal X (Thompson et al, 1997) with further 
adjustments by eyes and reference to two predicted sec- 
ondary structure models for SSU rRNA; the generalized 
eukaryotic SSU rRNA model (Peer et al, 2000) and 
the peloridiid 1858 RNA model (Ouvrard et al, 2000) . 
Some highly variable regions, which were not able to be 
aligned unambiguously across all taxa, even consider- 
ing the secondary structure, were excluded from further 
analyses. The aligned data are available from the au- 
thor or may be downloaded from the following web site: 
http: //life .zsu . edu . en/insect.185 7. 

Phylogenetic analyses were done by PAUP * 4.0b10 
(Swofford, 2002) and MEGA 2.0 (Kumar et al, 
2001). The method of constructing trees were NJ, MP, 
and ML. Unweighted parsimony analyses of various 
datasets were performed without considering gaps. As 
the number of taxa and the size of the data matrix often 
precluded more thorough searches, heuristic searches 
were performed with 100 random-taxon-addition repli- 
cates, and TBR (tree bisection-reconnection) branch 
swapping. When maximum parsimony tree was not the 
only tree used, either the strict consensus tree or 50% 
majority-rule consensus tree would be calculated. Clade 
stability was evaluated using two different parameters: 
bootstrap (Felsenstein, 1985) and branch support 
(Bremer, 1994). Bootstrap values were generated in 
PAUP " from 1 000 replicates, each with ten random- 
addition heuristic searches. Branch support values (a. 
k.a. Decay indices) were estimated with the program 
Autodecay V5.0 (Eriksson, 2001). We used the com- 
puter program MrModeltest 2.0 (Nylander, 2004) to 
identify the most appropriate substitution model for the 
ML analysis. Based on the result from the program, the 
ML tree was reconstructed using PAUP” . 


2 Results 


2.1 Sequences variation 

The length of the collected fragments ranged from 
942 base pairs (bp) in Maevius indecorus to 1 925 bp 
in Megacopta cribraria. The aligned region used in 
analyses comprised 1 910 bp, of which 793 (41.52% ) 
sites were variable and 519 (21.17%) sites were parsi- 
mony informative. The mean nucleotide composition 
was T 23.8%, C 23.8%, A 25.1%, G 27.3%. No 
significant bias in base composition among taxa were 


found. 
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Current family classification 


Tab. 1 Taxa used in this study 


Taxa 


Accession No. (GenBank) 


Length (base pairs) 





Aradoidea Aradidae 


Pentatomoidea 


Scutelleridae 


Tessaratomidae 


Plataspidae 


Pentatomidae 


Acanthosomatidae 
Canopidae 
Cydnidae 
Dinidoridae 


Urostylidae 


Coreoidea Coreidae 


Alydidae 


Hyocephalidae 
Rhopalidae 


Lygaeoidea Lygaeidae 


Rhyparochromidae 


Berytidae 
Piesmatidae 


Pyrrhocoroidea Largidae 


Pyrrhocoridae 


Idiostoloidea 


Cimicomorpha; Reduviidae 


Miridae 


Mezira granulata 


Mezira sayi 


Dysodius lunatus 

Coleotichus costatus 
Poecilocoris latus (Dallas) 
Cantao ocellatus (Thunberg) 
Tessaratoma papillosa (Drury) 
Megacopta cribraria (Fabricius) 
Aponsila montana (Distant) 
Erthesina fullo (Thunberg) 
Austrotechus rugosus 
Graphosoma lineatum 
Rhaphigaster nebulosa 
Elasmostethus sp. WCW-2003 
Stauralia compuncta 

Canopus sp. WCW-2003 
Allocoris sp. WCW-2003 
Megymenum sp. WCW-2003 


Urostylis westwoodi 


Notobitus meleagris (Fabricius) 
Mictis fuscipes (Hsiao) 

Cletus punctiger (Dallas) 
Anoplocnemis binorara (Distant) 


Cletus trigonus (Thunberg) 


Homoeocerus unipunctatus (Thunberg) 


Cletus sp. WCW-2003 
Leptoglossus occidentalis 
Leptocorisa sp. 

Leptocorisa acuta (Thunberg) 
Riptortus linearis (Fabricius) 
Riptortus sp. WCW-2003 
Maevius indecorus 

Harmostes sp. WCW-2003 
Serinthea sp. WCW-2003a 
Spilostethus hospes (Fabricius) 
Spilostethus hospes (Fabricius) 
Nysius sp. (Dallas) 
Graptostethus sp. 

Neacoryphus sp. WCW-2003 


Pseudopachybrachius guttus (Dallas) 


Paraeucosmetus pallicornis (Dallas) 


Paromius piratodes (Costa) 
Udeocoris nigroaeneus 
Neoneides muticus 

Jalysus spinosus 

Mcateella sp. WCW-2003a 
Mcateella sp. WCW-2003b 
Largus sp. WCW-2003 


Physopelta quadriguttata (Bergroth) 


Physopelta gutta (Burmeister) 
Dysdercus sp. WCW-2003a 


Dysdercus poecilus (Herrich-Schaeffer ) 


Trisecus pictus 
Emesaya brevipennis 


Lygus hesperus 


Campyloneura virgula 


AY252221* 
AY252222* 


AY252138* 
AY252274* 
AY627311 
AY627316 
AY627312 
AY627313 
AY627314 
AY627315 
AY252273* 
U88339° 
X89495% 
AY252322° 
AY252269* 
AY252269* 
AY252269* 
AY252224* 
AY252207* 


AY627321 

AY627317 

AY627323 

AY627329 

AY742883 

AY742886 

AY252261* 
AY252225* 
AY627320 

AY627320 

AY627320 

AY252161* 
AY252214* 
AY252226* 
AY252162* 
AY627319 

AY627325 

AY627326 

AY742885 

AY252411* 
AY627327 

AY627328 

AY742881 

AY252262* 
AY252412* 
AY252125* 
AY252164* 
AY252260* 
AY252227* 
AY742882 

AY742884 

AY252197* 
AY627318 

AY252213 

AY252321* 
U06476^ 


AY252317* 


1 


fa ee ey 
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— 





1 
1 
1 
1 
1 
1 
1 
1 
1 
1 


1 


021 
991 


058 
879 
866 
838 
923 
925 
924 
835 
883 
886 
924 
859 
877 
885 
863 
865 
104 


774 
804 
802 
813 
848 
859 
094 
875 
865 
865 
865 
016 
942 


1 889 
1 014 


831 


1 865 
1 823 
1 848 
1 894 
1 844 


988 


1 857 
1 879 
1 890 
1 048 


010 
101 
502 
835 
846 
100 
829 
038 
865 
922 


880 


Classification follows Henry (1997) .* Direct submission to NCBI GenBank database by Wheeler & Schuh (2004) b Campbell et al, 1994. *Aleshin 
et al, 1995. ^ Direct submission to NCBI GenBank database by Chalwatzis & Zimmermann (1996) . 
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Distance values among superfamilies and within 
superfamilies were determined by discarding positions 


with gaps in pairwise comparisons using the Kimura 


(1980) two- parameter model (Tab. 3). 


Tab. 2 Primers for amplification and sequencing 





Primer name Direction Primer sequence 
F Forward 5'-TCCCTGGTTGATCCTGCCAGTA-3’ 
R Reverse 5'-TAATGATCCTTCCGCAGG TTCA-3' 
2F Forward 5'-GGGAGGTAGTGACAAAAAATAACG-3' 
2R Reverse 5'-CCTGTTATTGCTCAATCTCGTG-3’ 
3F Forward 5'-GGTGAAATTCTTGCA TCGTC-3' 
3R Reverse 5'-ACATACTTGGCAAATGCTTTCGC-3" 
4R Reverse 5'-GTTAGAACTAGGGCGCTA TCTG-3' 


Tab. 3 Sequence divergence using Kimura two-parameter model 


Pentatomoidea 


Coreoidea 


Pyrrhocoroidea Lygaeoidea 





Variable range in 
a 0.4-3.4 
superfamilies ( % ) 


0.1-4.5 


0.7-7.7 0.4- 10.9 





Mean value between in-group species 4.6% 


2.2 Phylogenetic relationships 

The phylogeny of Pentatomomorpha was recon- 
structed using 185 rDNA. Trees recovered by MP, ML 
and NJ methods were mainly congruent (Fig. 1, 2, and 
3). 

Unweighted parsimony analyses yielded 24 equally 
parsimonious trees of 1 830 steps (CI 2 0.560, RI = 
0.624) ; the strict consensus tree is presented in Fig. 
1. The monophyly of infraorder Pentatomomorpha was 
strongly supported. The monophyletic Aradoidea ( Mezi- 
ra granulata, Mezira sayi and Dysodius lunatus ) was 
recovered by three methods. It was a sister group to the 
remainder of the Pentatomomorpha ( Trichophora ) , 
which was at a basal position of the tree topology. The 
monophyletic Trichophora was subgrouped into two ma- 
jor lineages. One of the two major lineages uncovered 
by 18S rDNA analysis comprised superfamily Pentato- 
moidea and Pyrrhocoroidea. Pentatomoidea was found 
to be monophyletic with strong support, which formed a 
sister relationship to Pyrrhocoroidea, except in the NJ 
tree (Pyrrhocoroidea being polyphyletic) . However, the 
family relationships of Pentatomoidea in three tree 
topologies were incongruent. The other major lineage 
recovered comprised Lygaeoidea, Coreoidea, and Idios- 
toloidea. The superfamilies Lygaeoidea and Coreoidea 
were both polyphyletic. In Lygaeoidea, four families 
(Piesmatidae, Berytidae, Lygaeidae and Rhyparochro- 
midae) were analyzed, grouping into two clades. One 
clade was Lygaeidae and Rhyparochromidae, which re- 
covered as a sister group with strong bootstrap value. 
However, the two families, Piesmatidae and Berytidae 


Variable range between in-group 


Variable range between 


species 0.396 — 12.7% superfamilies 0.3% — 7.8% 


appeared to have a sister relationship and were at a 
basal position of Trichophora lineages, distant to the 
clade of Lygaeidae and Rhyparochromidae. The super- 
family Coreoidea received low branch support in all 
trees, and is polyphyletic. The superfamily Idios- 
toloidea, represented by one spcies Trisecus pictus , 
seemed more related to the family Rhopalidae (super- 
family Coreoidea) . 


3 Discussion 


3.1 18S rDNA marker 

Our sequence based on phylogenetic analysis con- 
firmed the monophyly of Pentatomomorpha, Aradoidea, 
and Pentatomoidea. This was congruent with morpho- 
logical studies (Schuh, 1979; Wheeler et al, 1993; 
Henry, 1997) and the research of Li et al (2005). The 
Pyrrhocoroidea was further divided from the Coreoidea 
(s. lat) in this research, concordant with the results of 
Xie et al (2005) . Therefore we suggest that 18S rDNA 
was a proper marker to reconstruct the phylogeny of 
Pentatomomorpha. The placement of the Aradoidea 
within Pentatomomorpha was strongly supported by ML, 
MP, and NJ analysis, which was the sister group to Tri- 
chophora. In morphological studies, this is supported 
by the presence of pulvilli in some members (Tullgren, 
1918), tubular salivary glands (Southwood, 1955), 
egg structure, and trichophoran-type spermathecae 
(Pendergrast, 1957). The important implication of our 
analysis with respect to the higher taxonomy of Pentato- 
momorpha, was that Pentatomoidea was monophyletic 
and thus in accord with cladistic analysis of morpholog- 
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Lygus hesperus 
Campyloneura virgula 
Emesaya brevipennis 
Dysodius lunatus 

Mezira sayi 

M. granulata 

Mcateella sp. WCW-2003a 
M. sp. WCW-2003b 
Neoneides muticus 

Jalysus spinosus 
Harmostes sp. WCW-2003 
Serinthea sp. WCW-2003a 
Trisecus pictus 
Leptoglossus occidentalis 
Maevius indecorus 
Leptocorisa sp. 

L. acuta 

Riptortus linearis 

R. sp. WCW-2003 
Homoeocerus unipunctatus 
Cletus sp. WCW-2003 

C. punctiger 

C. trigonus 

Anoplocnemis binorara 
Mictis fuscipes 

Notobitus meleagris 
Udeocoris nigvoaeneus 
Paraeucosmetus pallicornis 
Pseudopachybrachius guttus 
Paromius piratodes 

Nysius sp. 

Graptostethus sp. 
Neacoryphus sp. WCW-2003 
Spilostethus hospes 
Graptostethus quadrisignatus 
Physopelta quadriguttata 

P. gutta 

Largus sp. WCW-2003 
Dysdercus sp. WCW-2003a 
D. poecilus 

Megymenum sp. WCW-2003 
Tessaratoma papillosa 
Allocoris sp. WCW-2003 
Coleotichus costatus 
Poecilocoris latus 

Cantao ocellatus 

Canopus sp. WCW-2003 


Austrotechus rugosus 


Elasmostethus sp. WCW-2003 


Stauralia compuncta 
Megacopta cribraria 
Aponsila montana 
Urostylis westwoodi 
Erthesina fullo 
Rhaphigaster nebulosa 
Graphosoma lineatum 


of Pentatomomorpha with 21 families and 53 species 
Bootstrap values are shown above nodes. Decay indices are shown below nodes. Current superfamily and family taxa following Henry (1997) are 


indicated on the right . 


Cimicomorpha 


Aradidae Aradoidea 


Piesmatidae 





Lygaeoidea 
Berytidae 





| Rhopalidae Coreoidea 


| Idiosroloidea 
| Coreidae 


| Hyocephalidae 


Alydidae 


Coreoidea 


Coreidae 


Lygaeoidea 


Lygaeidae 


Largidae 
Pyrrhocoroidea 


Pyrrhocoridae 


| Rhyparochromidae 


| Dinidoridae 

| Tessaratomidae 

| Cydnidae 

| Pentatomidae (Part) 


Scutelleridae 





| Canopidae 
| Pentatomidae Pentatomoidea 


Acanthosomatidae 





Plataspidae 





| Urostylidae 


Pentatomidae (Part) 


The phylogenetic trees using the Minimum parsimony (MP) method according to 18S rDNA sequences 
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Lygus hesperus 
Campyloneura virgula Cimicomorpha 


Emesaya brevipennis 














Dysodius lunatus 
Mezira sayi Aradidae Aradoidea 
M. granulata 
Mcateella sp. WCW-2003a |n T 
iesmatidae 
M. sp. WCW-2003b Lvelbsi dcs 
Neoneides muticus . 
. Berytidae 
Jalysus spinosus 
Harmostes sp. WCW-2003 Coreoidea 
Rhopalid 
Serinthea sp. WCW-2003a eae 
Trisecus pictus | Idiosroloidea 


Homoeocerus unipunctatus 

Cletus sp. WCW-2003 

C. punctiger Coreidae (Part) 
C. trigonus 

Leptoglossus occidentalis 


Maevius indecorus | Hyocephalidae Coreoidea 
Anoplocnemis binorara 
Mictis fuscipes Coreidae (Part) 
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ical data (Henry, 1997). 


phyletic Pentatomoidea differs from the three superfam- 


Morphologically, mono- 


ilies Pyrrhocoroidea, Lygaeoidea, and Coreoidea of 
Trichophora in many ways. These include the number 
of trichobothria, fusion and loss of parts of the plate- 
like ovipositor, an extra antennal segment, and a char- 
acteristic swelling of the spermathecal duct (Schaefer, 
1993). However, based on the tree topologies the rela- 
tionships among families of Pentatomoidea were much 
less ordered, which was consistent with the morpholog- 
ical studies (Stys & Kerzhner, 1975). Thus 18S rD- 
NA did not provide sufficient information to separate 
families inside this superfamily, suggesting requirement 
of more sequence data in phylogentic studies. 
3.2 Position of Piesmatidae and Berytidae 
Piesmatidae is firmly placed in the Pentatomomor- 
pha based on male genitalia, abdominal trichobothria , 
and pretarsal pulvilli (Drake & Davis, 1958), however 
its position in the suborder is inconsistent. It has been 
assigned as one superfamily, Piesmatoidea of Pentato- 
momorpha ( Sys, 1964, 1967; Stys & Kerzhner, 
1975; Schuh, 1986; Henry & Froeschner, 1988; Li 
& Zheng, 1994). Piesmatids have generally been 
placed in the Lygaeoidea (Schaefer, 1972, 1975, 
1993; Schuh & Slater, 1995; Henry, 1997; Xie et 
al, 2005). Henry (1997) grouped the Piesmatidae to 
superfamily Lygaeoidea by carinate pronotum, punctate 
hemelytra, the loss of inner laterotergites, and the bar- 
rel-shaped antennal segment | . However, these char- 
acters are not synapomorphous. Firstly, the barrel- 
shaped antennal segment I is also present in 
Aradoidea and in some species of Lygaeoidea. Second- 
ly, Aradoidea, Pentatomoidea, and many species of 
Coreoidea and Lygaeoidea have carinate pronotum. 
Thirdly the loss of inner laterotergites is quite homopla- 
sious, indicating a more detailed interpretation of their 
evolution (Henry, 1997) . The Piesmatidae are unique 
in having uniformly punctate hemelytral areoles, and 
two tarsomeres on each tarsus (also in Aradoidea and 
family Acanthosomatidae) while Lygaeoidea have three- 
segmented tarsi. Most significant is the presence of an 
m-chromosome (Jacobs, 1989) found only in the 
Coreoidea and most Lygaeoidea, however lost in a few 
taxa, e. g. Piesmatidae and Berytidae. Piesmatidae 
only has one pair of trichobothria on the abdomen 5 - 
6, and has lost it on the abdomen 4 — 5 (the same to 


some species of Berytidae). Li et al (2005) failed to 
reject that Piesmatidae is a family of Lygaeoidea. How- 
ever, based on the combined analysis of 18S rDNA 
and COX1, the MP and ML tree showed that Piesmati- 
dae and Berytidae were distant to Lygaeidae and Rhy- 
parochromidae (Li et al, 2005). Combining our 
molecular analysis, Piesmatidae was distantly related to 
other members of the Lygaeoidea; and according to 
morphological analysis, we suggest that the piesmatids 
should be assigned as a superfamily of Pentatomomor- 
pha rather than a family in Lygaeoidea. Morphological- 
ly, Berytidae have been placed in Lygaeoidea, along 
with the Clolbathristidae, Lygaeidae, Malcidae, and 
sometimes Idiostolidae (Sys, 1961) and Piesmatidae 
(Schaefer, 1993; Henry, 1997). 
sults of this study showed that Berytidae never merged 


However, the re- 


with the other families of Lygaeoidea and had a sister 
relationship to Piesmatidae instead. The reconstruction 
of their relationship will require further studies with 
other genes and more taxa samples. 
3.3 Polyphyletic Coreoidea and Coreidae 

With regard to the phylogeny of Coreoidea, there 
are many arguments in the morphological research. 
Henry (1997) and Li (1996) described the mono- 
phyletic Coreoidea with morphological data, whereas 
our analysis and result, derived from Xie et al (2005) 
and Li et al (2005), failed to describe the Coreoidea 
as monophyletic. This reflects the incongruence be- 
tween morphology and molecular data. For the family 
Coreidae, which is a large worldwide family, tradition- 
al classification is primarily based on the concepts of 
* synthetic systematics” (Hsiao et al, 1977; Schuh & 
Slater, 1995). The lack of morphological synapomor- 
phies suggests that this family is polyphyletic. Al- 
though without strong support, the tree topologies of 
our results also rendered this conclusion. The molecu- 
lar data and morphological data both indicate that tax- 


onomy of this family should be revised subsequently . 
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